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Abstract. The long-period perturbations in the oroit of the Lageos satellite due to the earth’s 
albedo have been found using a new analytical formalism. The earth is assumed to be a sphere 
whose surface diffusely reflects sunlight according to Lambert’s law. Specular reflection is not 
considered. The formalism is based on spherical harmonics; it produces equations which hold 
regardless of whether the terminator is seen by the satellite or not. Specializing to the case of a 
realistic zonal albedo shows that Lageos’ orbital semimajor axis changes periodically by only a few 
millimeters and the eccentricity by one part in 10-\ The longitude of the node increases secularly 
by about 3 x 10 - ^ arc sec yr — The effect considered here can explain neither the secular decay 
of 1.1 mm d <y - l in the semimajor axis nor the observed along-track variations in acceleration of 
order 2 x 10 - ^ ms'^. 



1. Introduction 

This paper is concerned with an analytical treatment of what is commonly called “the albedo prob- 
lem”, and how it relates to the orbit of the Lageos satellite. 

The basic idea behind the albedo problem is this: sunlight is diffusely reflected by the earth up 
to an orbiting satellite. The radiation pressure of the incident light exerts a force on the satellite. 
The task is to find this force and its effect on the satellite’s orbit. For an early review of the albedo 
problem, see Smith (1970). For some specific previous papers see Lochry (1966), and those cited 
in the next few paragraphs, as well as the ones cited by Smith (1970). 

The albedo problem is important because reflected sur.'ight significantly perturbs the orbits of 
satellites with high surface-to-mass ratios, such as the balloon satellite Page os (Prior, 1970; Laut- 
man, 1977a, 1977b; and Zerbini (1980). 

It is also important to know how reflected sunlight affects the Lageos satellite, even though it is 
small and dense. A high-precision orbit is necessary for Lageos’ mission of accurately measuring 
tectonic plate motion, polar motion, and length-of-day (Snith and Dunn, 1980). Smith (1983) has 
even speculated that the unexplained variations in along-track acceleration which have magnitude 
~2 x 10 -12 ms -2 (e.g., Rubincam and Weiss, 1985) are due to reflected sunlight. Anselmo et a!. 
(1983) support this view, while Walch (1982) finds this mechanism to be about an order-of-magni- 
tude too small. Given this state of affairs, it is worthwhile to look further at the albedo problem 
and Lageos. 

Solving the albedo problem analytically is not easy. Even though irradiance falls off with the 
inverse square of the distance, like gravity, the albedo problem has a number of complications not 
found in the gravitational problem: (a) The irradiance on the satellite depends on the orientation 
of the surface element doing the reflecting as well as the distance; an example is Lambert’s law 
(e.g., Brown, 1965, p. 225); (b) The irradiance also depends on the solar zenith angle; (c) Only sur- 
face elements illuminated by the sun and visible from the satellite contribute to the force; (d) The 
earth’s albedo varies with time, as well as latitude and longitude; (e) The atmosphere reflects and 
scatters light, so that one is not dealing with just reflection off the earth’s surface; (0 The force 
depends on the properties and orientation of the satellite; general!'’ a fiat mirror will feel a dif- 
ferent force from a mottled ellipsoid, for example. 

Given these complications, it is not surprising that Zerbini (1980) and Walch (1982) forsake the 
analytical approach of the other papers mentioned above for a numerical one. However, there are 
two good reasons for pursuing an analytical treatment: (i) equations give more insight into the 
problem than numerical methods, and (ii) equations can often lead to significant savings in com- 
puter time compared to numerical methods when integrating satellite orbits over time. Thus the 
rationale for the present work. 
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The complications (a) - (0 mentioned above are dealt with here as follows. The earth is as- 
sumed to be a sphere whose surface diffusely reflects light according to Lambert’s law. The al- 
bedo variation on this surface is assumed to have a general spherical harmonic expansion whose co- 
efficients may be functions of time. Spherical harmonics are also used to deal with solar illumina- 
tion and the spherical cap seen by the satellite. The satellite itself is assumed to be a sphere of uni- 
form albedo. 

These assumptions lead to a vector equation for the albedo acceleration of the satellite. Then 
expressions for the orthogonal accelerations R. S, and W are found after specializing to the case of 
a zonal albedo variation (i.e., the earth albedo depends only on latitude). These accelerations are 
to be used in the Gaussian form of Lagrange’s planetary equations; the acceleration equations have 
the virtue that they hold regardless of whether the satellite sees the terminator or not. Specializing 
still further to the case of Lageos and a realistic zonal albedo model based on Stephens et al. (1981) 
yields the following long-period results. The semimajor axis changes periodically with an amplitude 
of ~5 mm. The eccentricity changes periodically by about one part in 10^. The longitude of the 
node increases secularly at the rate of about 3 x 1G - ^ arc sec yr - *; periodic changes have amplitude 
~ 3 x l0 -< * arc sec, as do the changes in inclination. All of these effects are quite small. Lambert- 
ian reflection cannot account for the secular decrease of 1.1mm day - * in the semimajor axis or the 
observed variations in along-track acceleration which have amplitude ~2 x 10 - *^ ms - 

Some notation is given in Tables I and II; other quantitities are defined in the text. A photo- 
graph of Lageos is shown in Figure 1. 

2. Derivation of equations 

In the case of the gravitational field, the usual method of obtaining the change in the Keplerian 
elements of the orbit with time is to first find the disturbing function and then use it in Lagrange’s 
planetary equations (e.g., Kaula, 1966, pp. 25 - 37;Caputo, 1967, pp. 140 - 158). Unfortunately 
this approach cannot be used in the present problem where the reflected sunlight follows Lambert’s 
law. The reason: the force cannot be written as the gradient of a potential, as it can for the case of 
the gravitational field. This is proved in Appendix 1. Hence no disturbing function exists. Instead 
tiic force must be dealt with directly. Once the force is found, it can be used in the Gaussian form 
of Lagrange’s planetary equations, as is outlined below. 

The derivation of the magnitude and direction of the force on a satellite due to diffusely reflect- 
ed sunlight will now be given. Consider the inertial reference frame as shown in Figure 2. The 
origin is at the center of the earth. The frame does not rotate around with the earth. Suppose the 
sun shines on a surface element dS as shown in the figure. The sun is at position r<j while dS is 
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at r^. The surface element diffusely reflects some of the sunlight up to a satellite along r^ - r^, 
wheie tl is the satellite position. This light exerts an infinitesimal force dF on the satellite due to 

all of the sunlit surface elements visible to the satellite. This net force will be called the albedo 
force for short. Likewise the albedo force divided by the satellite mass will be called the albedo 
acceleration. 

The albedo force depends on the shape, reflectivity and orientation of the satellite as well as on 

the characteristics of sunlight and the albedo of the earth. However, if the satellite is a perfect 

sphere with uniform reflectivity, the force dF will always act along the surface eleme.it-satellite line, 

-*■ -*■ -*• 

in the direction away from dS; in other words, along r^ - r^ as shown in Figure 2. Such a case 
obviously greatly simplifies the treatment of the albedo force. Fortunately, Lageos satisfies these 
conditions to a very high degree of approximation (see Figure 1). Hence in the following the 
infinitesimal force dF will always be assumed in the dirction of r^ - r^. 

Now that the direction of dF has been established, its magnitude will be computed. Since the 
sunlight is assumed to be diffusely reflected according to Lambert’s law, the radiance of dS, mea- 
sured in watts per steradian per meter squared, is Nrcos/3 (e.g., Brown, 1965, pp. 220-225), where 
P is the colatitude in the local reference frame (Xj, y j, Zj) of dS (see Figure 3) and Nr is a constant 
to be found next. 

The radiant flux incident on dS is Fg cos dS, where Fg is the “solar constant” and tf/g is the 
colatitude of the sun in the local frame of Figure 3. The amount of flux reflected by dS is A(0, X) 
Fg cos \p g dS, where A(0, X) is the albedo of dS, and (0,X) are the colatitude and longitude, respec- 
tively, in the inertial frame of Figure 2. This flux must be equal to Nr cos 0 dS, integrated over the 
hemisphere into which dS diffusely reflects sunlight, assuming the power into the surface element 
equals the power out. Therefore 


2 n n/2 

dS / / N r cos P sin/J d0 da = A(0,X) Fg cos^g dS , 


o o 


where a is the longitude in the local frame of Figure 3, so that 
Nr = A(0,X) Fg cos tf/g/ff 
The power E incident on the satellite due to dS is 


E = NRCOS^dSdA = A(0,X) Fg cos ^g cos dA/ff, 

— ► — ¥ 

where ^ l the angle between r ^ - r^ and r^ (see Figure 2) and dA is the solid angle subtended by 
the satellite from dS. The magnitude of dfc can be found from E in the following manner. The inci- 
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dent momentum p of the light is related to its energy E by E = pc, where c is the speed of light (e.g., 
Alonso and Finn, 1968, p. 14). The incident momentum per unit time is thus E/c. The force on 
the satellite is the change in its momentum per unit time, and is given by Cr E/c (e.g., Rubincam 
and Weiss, 1985), where Cr is the radiation coefficient. (If the satellite absorbs all of the incident 
light and re-radiates it isotropically, then Cr = 1. If the satellite is a flat, perfect mirror reflecting 
the incident light straight back, then Cr = 2. So 1 <Cr<2. For Lageos, Cr = 1 .13). The magni- 
tude of dF is ihe'?fore 

- Cr F S 

|dF| = A(9,X) cos cos dS dA. (1) 


The vector force dF can be found by multiplying ( 1 ) by the unit vector (r^ - r A^l r L “ r Al - Also, 
for a spherical satellite of radius R^, dA it R^ / |r^ - r^ I to an incredible degree of accuracy, 

—4 -4 

since R^ « I r^ - r^| always. Further, 


cos * s = ( r S - r A' * r A/(lrs-^| • ^/(Irgl • lr^|) 

because |r^| « |r<j|, since the sun is so far away from the >’arth. Likewise, Fg * F§ (a§/|r§|)^, 

O 

where a^ is the semimajor ?xis of the earth’s orbit ^bout the sun and Fg is a constant, equal to 
about 1376 W m~ 3 (Hickey e_t al„ 1980). No such approximation can be used for cos \p so it is 

(fL ~ i^) * ix/dfL - f aI * |r A () exact| y- 

-4 

The albedo force F can be found by substituting the formulae above into (1 ) and integrating over 
area yielding 


?= M L r L = 


C R r L 2 f S° 


f 7 


A(0 ,X) [r s • r A l i(rL-rA> ,r AH r L- r A) 
I^L " ^A |4 |r A' 2 


dS (2) 


in vector notation where r L is the albedo acceleration in the reference frame of Figure 2, is the 
satellite mass, and Z is the sunlit portion of the earth visible to the satellite. The albedo accelera- 
tion can be written 
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r L = 


c r R L 2 Fs° 

cM l 



//, 


A(0,X) cos (cos 17 -p ) 


Z (1 - 2p costj + p 2)2 


• [ (sin 0 l cos X^ - P sinfl cosX )x 

+ (sin sin ^L - p s * n 0 s ‘ n ^ ) y 
+ (cos 0 l - P cosfl ) z] sin0 dfl dX 


( 3 ) 


after reverting back to trigonometric notation. Herer? is the angle betwen rj^ and r A (see Figure 
2), so that cos r? = r^ • r^ /(|r^| • |r^l),R A = ir^l, r L = |r L l, r s = |r s l, and P = Ra^L* A1so ’ 
uS = R a si.i 6 d0 dX, and (0^, X^) are the colatitude and longitude of the satellite respectively in 
the frame of Figure 2, whereas (x, y, z) are the unit vectors along the appropriate axes in the frame 
of Figure 2. A spherical earth of radius R A is assumed in (3); R A is the distance from the center of 
the earth to the “top of the atmosphere”. The “top of the atmosphere” is taken to be 10 km 
above the surface (the height of the troposphere) so that R A = 6371 km + 10 km = 6381 km, where 
637 1 km is the radius of a sphere with the same volume as the earth. (A more common choice for the 
top of the atmosphere is 30 km; e.g., Bess and Smith, 1981.) 

The task ahead is to make the integral in (3) tractable. Spherical hai monies form a viable ap- 
proach, as shown below. 

Suppose a function D is defined so that 

n 

fcos^^ 


D =< 


IT 

—< 4/ s < * 


then D can be written as an infinite sum of Legendre polynomials 


D = 


OO 


I 


d^ Pl ( cos tf'S ) 


where (cos 4/ §) is the Legendre polynomial of degree L and the d^ are the appropriate coeffici- 
ents. (The subscript should not be confused with that used for the satellite.) The addition theorem 
for spherical harmonics 
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P L< cos *S> * T777 2 Yl'*(# s ,x s i y l j (0 .XI 


+L 


2L-M j=_! 

(e.g., Merzbacher, 1970, p. 189) can be used in the above equation to write 

co +l 


D(0 ,x, e s 


’ x s ) = 4n 2s Z (®s <x s ) y l 

i ^o 


where the 


Y L (0,\) = (-1) 


J+JJI 
2 


(2L+1XL- I J |)! 

[ 4 jr(L+ |J | )! j 


V/2 


F L,|J| ( cos0) e 


iJX 


( 4 ) 


are the spherical harmonics used in quantum mechanics (Merzbacher, 1970, p. 185), and (0§, X§) 
are the respective colatitude and longitude of the sun in the frame of Figure 2. Here Pl j j j(cos 0 ) 
is the associated Legendre polynomial of degree L and order |J| and i = V^T. The asterisk (*) 
denotes the complex conjugate and -L < J < + L. The quantum-mechanical spherical harmonics 
have normalization 



(0.X ) 



(0,X ) sin0 d0 dX 


= 6 


JJ’ 


6 ll* 


(Merzbacher, 1970, p. 186) where 6 jj> is the Kronecker delta. The relation between the quantum- 
mechanical spherical harmonics and those commonly used in geophysics is given in Appendix 2. 
Likewise a function B defined so that 


r 


B = \ 


can also be written 


cos rj - p 

, , _ 7,7 

I ! -2 p cos i? + p J A 


0 


0 < tj <Arc cos (p ) 


A rc cos (p ) < r? < t r 


B = ^ b g p y (cosr? ) 
G=0 
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so that 




4rr 


2 

«=0 




(0L’*L> Y ™ <M) 


(5) 


where the bg depend on p. Arc cos (p) is the angular radius of the cap visible to the satellite. Find- 
ing the dL and the bg will be postponed until after a general expression for the albedo acceleration 
is derived. 

Now D(0, X, 0g, Xg) is zero over the night hemisphere of the earth. Also. B(0, X, d X^) is 
zero outside the spherical cap visible to the satellite. Hence by substituting D(0, X, 0g, Xg) for 
cos tf/g and B(0, X, 0^, X^) for (cos r? - p)/| 1 - 2p cos r\ + p- ] - in (3) allows area of integra- 
tion to be changed from * to the entire sphere, since now the integrand is zero outside the 
region X . This enlargement in the area of integration greatly simplifies the evaluation of the 
integral, as will be seen shortly. 

The albedo A(0, X) is also assumed to have a spherical harmonic expansion: 



N=0 K=-N 

where tne coefficients A^ refer to the inertial reference frame of Figure 2. The albedo is de- 
fined over the entire earth, regardless of day or night. This should cause no conceptual diffi- 
culties; the values associated with the night side may be thought of as those the albedo would take 
if the sun could somehow be made to shine on the sunless hemisphere. Further, 0 < A(0, X) < 1. 
Now realizing that 


cos X = 


e +iX + e -iX 
2 


sin X = 


c +iX _ e -iX 
2i 


and using the properties 
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cos0 Y™ ( 0 , X) 


y (6-m+l)(6+m+l) 
(26+!)(26+3) 


v m io . /(8--m)(g +ml\ ,m 

^+1 ’ \/ (26- 1 )(26+l ) Y ®“ 1 


. . +i\ ..m /(6+m+l )(8+m+2) 

S,ne C Y 6 <•*> - "V (26+1X26+3) 


n.+ l 
6 + 1 


(0.X) 


+ 



( 6-m-l ) ( 6-m) 
(26 — 1 ) (26+ 1) 



(0.X) 


sin0 e 


-iX 


,m 


Yn (0.X) = 


/Tu 

V (2 


( 6-m+l)({-m+2) „m-l 


(26+1 )(26+3) 


‘ 6+1 


(0.X) 




(6+m- 1 )(6+m) m-1 

(26 — 1X26 + 1 ) c -> ( ’ * 


(e.g. Merzbacher, 1970. p. 187) allows one to write from (5) 


B(0, \0 L .X L ) (sin0^ cos 0^ -psin0 cosX) = 2 v 


2 

6=0 


+6 

J-t 


Ye™ < fl L- X L } 


r 

< - 

k. 



+ 



- P 


U\ i 
\26 + l/J V 


! (6+m-l ) (6+m) ,,m- 1 _ 

Y„ , ( 0 . X) 

( 26—1 X 2 6 + 1 ) 6-1 


- p 


(H\ 
\ 2 * + 1 /_ 




(6- m-1) (6-m) Y m+1 (0X) 

(26 — 1 )( 26 +1 ) 6-1 


(0.X) 
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after gathering coefficients of Y™* (0^- X^). d nterRS -i n ^y< there is no danger of the superscripts 
or subscripts of the Y™ * * (0,>.) exceeding their proper limits in the above equation; the radical in 
front of such terms is zero.) 


Using the equation above, similar equations for B(0 ,X ,0 j_,X l) (sin 0 1 _ sin X “ P s ' n ® sin X ) and 
B(0,X,f ^ ccs 0 L ~ p cos 0 as wel * as and ’ n and > n,e 8 ratin 8 ovcr the ^hole sphere 
yields 




9 






(or the albedo acceleration, where 
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•<> 



2n v 

// 

*'o •'o 


*K J m 

Y n ( 0 ,X)Y l (0,X) Yg (0 ,X) sin0 d0 dX 


= ( 1)K j' (2N-H)(2L-H)(2g-H) j‘^ 


I n l e 

\-K J m. 


N L C 


0 0 0 , 


(e.g., Rotenberg, et al., 1959, p. 5) is the simplification of the integration mentioned earlier. The 
two-row arrays are the 3-j symbols of quantum mechanics (e.g., Rotenberg et al., 1959, p. 2) 


c 

\ m l 


j 2 i 3 


= ( _, ) j r j 2 _m 3 


(■ 


Vi 


0j + j 2 -i 3 ) ! (i, - j 2 ♦ j 3 )! (-j, + j 2 ♦ j 3 )! 0j + m,)! (J ! - ",)! °2 r m 2 )! °2~ m 2 )! °3 4 m 3 )! 0 3 ~ m 3 )! 

0l + i 2 + j 3 +1)! 


(- 1) 1 


^ k!(l 1 +j 2 -j 3 -W!0| -mj-k)! 0 2 ♦ m 2 - k)! (i 3 - j 2 + m, + k)! (j 3 - j, - m 2 + k)! 


The sum in this equation is over all values of k so as to make the right side meaningful. The 3-j 
symbol is automatically zero unless the triangle inequalities 


•>1 + J2 “ h > 0 

i\ ~ h+h > 0 

-Jl + J2 + >3 > 0 

are satisfied (Rotenberg et al., 1959, p. 3). The 3-j symbols also have the following properties: 


I h h 

M 

= 0 

for m j + m 2 + m^ f 0 

( 8 ) 

\ m l m 2 

m 3 y 
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so that 


(- 1 ) 




Jl J 2 J 3 
rnj m-> m 3 — 1 i 


= (-D 


Jj J 2 


m | m 2 m 3- 1 


for example, and 


I j i h h' 

l -mj -m 2 -m 3 


Jl+b+h 
= (-l) 



(9) 


(10) 


(Rotenberg et al., 1959, pp. 2-3). These properties will be used shortly. 

The albedo acceleration r£ of the satellite will next be expressed in terms of the orbit parame- 
ters. 

First assume that the albedo variations are attached to the eaiih, i.e., the light and dark markings 
rotate around with it. In this case the coefficients of ( 6 ) can be written 




e -iK e 


(ID 


K 

where the are now the earth-fixed coefficients and 0 is Greenwich sidereal time. The 
right-handed earth-fixed frame has its x-axis intersecting the equator at the Greenwich meridian 
and its z-axis piercing the north pole. 

It can easily be shown from the standard development of the gravitational potential (e.g., Caputo, 
1967, pp. 140-145; beware of typographical errors) that 


Y f m (0,X) 


m+|m| 
(- 1 ) 2 


^ F C|m|p (,) j 


p =0 


e i(26 m |m| -1) I(£-2p)(w+f) + |m|fi) C-m even 

^ 6 m|m| ~ 1 * i(26 m | m | -1) [(C-2p)(w+ f) + |m|f2] 
i 


C-m odd 


( 12 ) 
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In this equation Fg m p(I) is the normalized inclination function 




(2-6 om )(28+l)(8-m)! 


(£ + m>! 


^£mp ^ 


(Rubincam, 1982, p. 364), while I is the inclination of the satellite orbit with respect to the earth’s 
equator, ft the position of the orbit node on the equator, u> the argument of perigee, anu f the true 
anomaly. 

Using ( 1 1 ), ( 1 2), and the succeeding equation, one can write 


,nr 


a *Wv y * (om 


= (-!) 


K+|K|+J + |J| +m + |m| r (2-6 oK ) “f 4 L 8 _ 

^7 X S F L|J|h (1 S )F C|mlP U) 

4r^ L(2-6 oJ )(2-6 om )J h=Q p=0 

+i {(1-26 J( j|) [(L-2h)(u s +f s ) + |J|ft s ] + (l-26 m|m |)[(8-2p)(w+f) + imlft] + K 6 | 


/ A N|K|1 —iO — 25 K|K|> A N|K| 2 


L - J even 
£ -m even 


•< 


(1 ~ 25 m|m| ) A Nl K 1 1 “ < 1 ~ 25 Kl Kl )( 1 ~ 2h ml d > A Nl Kl 2 
i 

(1 ~ 26 J |J l } A N|K|1 - ( ,_26 K|K| )(1_26 J|J1 ) A N|K|2 
i 


II: 


L - J even 
£ - m odd 


L - J odd 
' £ - m even 


V - 26 J|J|) (, - 26 m|m| ) A N|K|1 + » (, " 26 K|K1 )(1 - 26 J|.l| )n “ 26 m|m| ) A N|K|2 


]- v ' ; L - J odd 
£ - m odd 


(13) 

where A^ug j and A^u^ ■> are the earth-fixed albedo coefficients of the geophysical spherical har - 
monies; these harmonics are defined in Appendix 2. Also, 1$, ft§, ojj, and fij arc the orbit parame- 
ters for the sun. 
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Now realizing that any sum of the following sort can be written 


+N +L + fi 

2 2 2 C KJm 

K=-N J=— L m=-i 


+N +L +6 i 

2 2 2[<^k><'V (,+ w J 

K=0 J=0 m=0 


{ C KJm + C KJ-m + C K- Jm +C K-J-m + C -KJm + C -KJ-m + C -K-Jm + C -K-J-m j 


and that 


K -K 

(-1) = (-D 


6 oj 6 o-j 


as well as making use of (8) - (10) and (13) in (7) yi ids, after a great deal of labor, the albedo accel- 
eration expressed in terms of the satellite orbit parameters: 


"Rl c r f s 


cM, 



oo +N °o +L °° +C 


1 1 1 11 1 

N=0 K=0 L=0 J=0 C=0 m=0 


(-1 ) f 


(2-6 oK ) (2N+1 ) 


- i '/ 2 


L (2-6 0 j)(2— 6 om )(2L+l K2C+1 )J 




. (I+5 oK ,,1+, oJ ,,I+8 om*J 


+L +C 

* 2 2 ^Uh (I S ) ^8mp (,) 

h=0 p=0 
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a NK1 cos < V LJh + V Cmp - K e )+ A NK2 sin(v£ Jh + Vg lip - K e ) 
A NK i sin(V^ lh + V^ p - K 6 )-A NK2 cos(v£, h + V e ^ p “ K ® ) 
A NK1 sin(V^ Jh + Vg mp - K 9 ) - A NK2 eos(VL Jh + V^ mp - K © ) 
-A NK1 cos (Vu h + Vg mp - K 9 ) - A NK2 sin (V^j h + Vg mp - K © ) 




I 

II 

III 
IV 
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A NKl COs(V u h - V 8m P +K€> )_ ^NK2 sin ( V LJh- V C L mp + Ke > 
- A nk , sin (V^, h - V^p + K 0 ) - A NK2 cos(v[j h - Vg L mp + K e ) 

A NK | sin (Vlj H - \^ mp + K 0 ) + A NK2 cos(V^j h - v £ L m p + K e ) 
ANK 1 Cos(V^ h -V c L mp+ K e )-A NK2 sin(V^j h -Vg L mp +K 0 ) 


I 

II 

III 
IV 
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S j 

f A NK1 cos(Vuh + ^mp +K 0 ) “ A N K2 s in(Vu h + V Km P + K 6 

s 

J A NKI sin (V S uh + V 8 L mp + K e ) + A NK2 cos(V Uh + ^ mp +K e 

A N Kl si "< V Uh + v e L mp +K e > + A N K 2 cos < V Uh + V (: L mp +K 0 

S T 

\-A N Kl cos ( v Lh + V Wp + K 0 > + %2 sin % h ^; mp +Ke 


) 1 

) II 

) III 
) IV 



J(i- m-ixe- 


m) 


' b 8-l\ _J M 

1/ P ^2C+1/ 


N L 8-1 


_K J m+l 




N L 8-1 


UK J m-1 1 




N L 6+1 

UK J m-1, 


- AwceosOLr^^p-K e ) - A N K 2 "n% h *V 8mp -K 0 ) 

- * N KI»%„ , i |) - |< 0 > - A N K2“"( V f,K tV Ln- K 9 > 


Uh 8mp 

I NKI an %H tV e m p‘ K 0 )+ A N K2“*( V u h tV Lp- K 9 * 
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S L 

( A NK1 sin (V S uh - V g mp - K 0 ) - A NK2 cos (V uh - Vj mp - K 0 ) 

) A NK ,cos(Vu h -Vc mp -K 0 ) + A NK2 sin(Vu h -V emp -K0 ) 11 

-A NK 1 cos( v L)h _ v'jmp- K 0 ) - A NK2 sln <V Uh “ V 6mp ” K 6 * 

, A NK ,sin(V S LJh -V^ p -K0 » - V 2 cos V^- K 0 , W 
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“• S L — § 

“ A NK1 s« n ( v Uh - v g mp + K 9 ) “ A NK2 cos < V LJh " V £ mp + K © ) • 

- A nk , cos (Vjj h - + K 0 ) + A NK2 sin (V^j h - + K © ) 11 

+ A NK1 cos (V Uh - V Cmp + K e ) - A NK2 sin (V^ - Vg L mp + K © ) III 

•i j S L 

~ A NK1 sin ( V LJh~ V Cmp + K © ) - A NR2 cos (V LJh - V £mp + K © ) IV 
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- A NK1 sin(V Uh + V Cmp + K 9 ) - A NK2 cos(V Uh + V 2mp +K 9 ) 
r* L * S I 

+ A NK1 cos ^ V LJh + v emp + K 9 ) “ a NK 2 sin (V LJh + V Cmp + K0 ) 

A NK1 cos(V Uh + V Cmp + K 0 } " A NK2 sin(V Uh + V Cmp +K 0) 

A NKl sin(V LJh + V Cmp + K 6 > + A NK2 C0S < V Uh + V Cmp + K e) 




— S L - S L 

A NK1 cos(V Uh + V Cmp - K 9 ) + A NK2 sin (V LJh + V cmp “ K e) 

— S L - S L 

A NK1 sin(V LJh + V Cnip“ K 9 ) ~ A NK: cos(V LJh + V Cmp" K e) 

A NK1 sin ( V Uh + V 8mp “ K 9 ) “ A NK2 C0S (V LJh + V Cmp _K 9 ) 

- NK1 tos( v uh + ^mp' K 9 } - A NK2 sin(V LJh + V fimp _K 0 } 
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A NK1 cos(V LJh - V fimp ~ K 6 ) + A NK2 sin(V LJh V £mp 


-K 0 ) 


S L S L 

- A NK , sin(V Uh -Vj mp -K * ) + A NK 2 cos(V LJh - V {mp -Kel 

S L — S L 

a NK 1 sin ( V LJh ~ V £mp - K 0 ) - A NK2 cos(V LJh~ V 8mp- K 0 > 
S L ™ S L 

A NKl cos ( V Uh- V Cmp- K 0 ) + A NK2 sin ^ V LJh- V cmp- K 0 ) 


+(-l) 



£-1 


0 


1 or 

1 re 

1 

(-1 


/(£-m)(£ +m) 

\2£ -1/ 

\2£+l / 

V 

L' 




N L 
\- K - J 


I 

Ii 

III 

IV 

£ + 1 
m , 


N 


0 0 



- P 


(-) 

\ 2£ + l/ 


y(£-m+l)(C+m-H) 



NK1 


cos 


( v Uh- V Cmp + K 0 > - A NK2 sin(V Uh- V £mp + K0 ) 


- A 


NK1 


sin( v LJh- V £mp + K 0 > " A NK2 cos (V Uh ~ V £ mp + K 0 > 


S l » S L 

A NK1 sin(V Uh " V £mp + K 0 ) + A NK2 cos(V Uh ■ V £mp + K 0 ) 

S L — S L 

a NK 1 cos(V Uh “ V £mp +K 0 ) ~ A NK2 sin(V Uh " V £mp + K 0 ) 


I 

II 

III 

IV 


-H-l) J+m < 


N 


L 

0 



b £-l 

2£-l 



N 

-K 





^(E-m+l )(£ +m+l ) 
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A NK 1 c °s( v LJh + V 6 mp + K 0 > - A NK2 sin(V Uh + V fmp + K 0 > 1 

A NK1 sin < V LJh + V 8mp + K 0 > + A NK2 cos (V Uh + V £mp + K © > 11 

_ S L S L ' L >( I4 > 

A NKI sin(V LJh + V fmp + K 0 > + A NK2 cos(V Uh + V £mp + K Q > 1,1 I 

- a NK1 cos(V LJh + V «mp +K 0 } + A NK2 sin(V Uh + V 8mp + K 0 > IV J 1 


In this equation 


V LJh = (L - 2h)(^ s + f s ) + Jft$ 


V Cmp = * ®~2p)(cj + 0 + mft 


It should be noted that K, J, and m are now non-negative numbers, as may be seen from the summa- 
tion symbols in (14). Also, only the Roman numerals for the four conditions are given to save from 
writing "L - M even, C - m even ” etc., over am’ over. 

It is of some interest that the leading subscript on the inclination functions can take on the values 
0 and I , in contrast to the gravitational case. These functions are listed in Table III, since they are 
not given in the standard references. 

The albedo acceleration can be expressed in terms of the orthogonal components (R, S, W) used 
in Gauss’ form of Lagrange’s planetary equations. These components are derived via a rotation 
matrix from the Cartesian components, (x. y. x) found in ( 14): 



cos (cj + 0 cos ft 
-cos I sin ft sin (cj + 0 

-sin (u> + f ) cos ft 
-cos I sin ft cos (u> + f) 


sin 1 sin ft 


cos (cj + D sin ft 
+cos I cos ft sin (cj + 0 

-sin (cj + f) sin ft 
+cos I cos ft cos (cj + f) 

sin I cosft 




(15) 


(Rubincam, 1982, p. 370). In this equation R is the radial component of acceleration, S is the trans- 
verse component in the plane of the orbit, perpendicular to R and positive in the direction of motion, 
while W is normal to the orbital plane and positive in the direction of the angular momentum (e.g., 
Blanco and McCuskey, 1961, pp. 177-178). Of course r^ = "x* x + y* y +7 z. 
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The general equations for R, S, and W will not be written here. Instead only the case of zonal 
albedo variations will be considered, so that the light and dark markings of the earth vary only with 
latitude. The reason for doing this is that only the zonals will give a long-term change in the Kep- 
lerian elements, barring a resonance. So K = 0 from this point onward. 

Using the general relations 


cos a cos b 
sin a sin b 
sin a cos b 
cos a sin b 


Vicos (a+b) + Vicos (a -b) 
-Vicos (a+b) + Vicos (a-b) 
Visin (a+b) + Visin (a-b) 
'/isin (a+b) - Visin (a-b) 


(16) 


(e.g., Kaula, 1966, p.31) together with (14) and (15) yields the following equation for the radial 
albedo acceleration: 


R = 


ffR L C R F S 


4cM L 



OO UO -+L oo +8 


2 2 2 2 2 

N=0 L=0 J=0 6=0 m=0 


d L A NOI 


(2N+I ) 


Vi 


(l-(-l) N+L+C l 2 +L +6 _ 

2 2 

h=0 p=0 


L(2-6 oJ K2-6 om )(2L+l)(2B+l)J ^ ' Uh^S* F «mp (I) 


/( 8-m+l )(8-m+2) (— ] -pi— \ 

V L\^ + 3/ \28+l/ 


-I /N L 8 + l\ / N L 8+1 ^ 


- /( g+m— 1 )(K + mi 

V 


f D 6-l 

128-1 


(26 + 1 


^0 0 0 / \0 
r N L 8 — 1\ / N 


-cos (Vy h + V j mp - (w+fHSl )] 
. <! -sinl v | J h + vl {mp-^+Wll | 
-»n ( V^jj, + v)r m p - 
♦cos [ vfj h + V t L mp - )| 


J m- 1, 

l e-i 1 


(0 0 0 / \0 J m-li 

I 

II 

III 
IV 


1 +COS I 
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-COS I V S + V S + (Wirf -n )] \ 

LJh £mp 

-sin (V s , + V S + (w+f- ft )] 11 

Uh C mp 

-sin (V^jh + V^p + M-f -ft )] 111 

+ cos(V^j h + vJ mp + (uHf-n)] IV 



+ co.lv^ h + v £ L mp — ( W ff-n)| i 

rsin l V^jh * V^mp “ (tJ+f- ^ ) l 11 

+sin I v Uh + V Cm P -<w +f -n)l 111 

-c°s|Vuh + v }'mp- <u,+f - n )l IV 
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( +cos[Vu h + V £ I ; np +(u;+f-ffl)] 1 

1 + sin l V lj h + V g L mp + k+f+fi)] II 

\ S L 

I sin (Vi iu + V +(j+f+J2)) III 

I LJn Emp 

( -cos [vjj h + V^+ (w+f+n)l IV 



-cos [V^, h - Vg L mp + (w+f-S2)l I 

+ sin (V^j h - V c L mp + (orhf-n)] II 

-sin [V^jh - Hirr + M'f -0 )1 HI 

-cos [V^j h - V^p + fcj+f-12)] IV 



-cos (V^ - \^ mp - (ort-f+ o)) I 

+ sin (V?i h - V L - (w+f+n )] II 

LJn Emp 

- sin ( v? Jh - v c L mp - M* )! HI 

- cos ^uh-^p-^f+nH ,V 
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+cos,v uir v « L mP + ( " +f+n)l 

S L 

-sin l V , " V 0 + 

LJh fimp 

S L 

+sin (V -V + (w+f + ft)] 

LJh Cmp 

S L 

+cos [ V - V + M- i )) 

LJh Cmp 


I 

II 

III 

IV 



, +cos | V 
-sin IV 
•fsin (V 
+ cos [V 


S 

LJh 

S 

LJh 

S 

LJh 

S 

LJh 




<w+ r-nii 

(uj+f- n )i 
(CJ+ f- 11)1 

(w+t n)i 


i 

ii 

in 

IV 
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+ sin [VLJ h + v 8mp + (u;+ ™ 

~ cos l v Uh + v 8mp + ( w+ 01 
S L 

—cos [V + V +(u+01 
LJh C mp 

-sin [Vuh + V C mp + (o;+f)1 


I 

II 

III 

IV 





S L 

“ sin (Vy^ - V^p 

S L 

- c°s( v uh " V 8mp 

S L 

+ cos[Vjj h - Vg mp 

S L 

-sin [V LJh - V Cmp 


- (u) + 01 


- (w + 01 


(Oi + 0) 


- (u) + 0) 


+2(— 1 ) m N C-m+1 )(C +m+l ) 


N L e + l\ / N L CH 


0 0 0 / \0 J -m 


+ V( c -m)(C+m) - p 

I \ 2C — 1 / \2fi + l 


N L «-l\/N L C-l 


0 0 0 / VO J -m 


4-sin l v LJh- V g mp + (w+ 01 

+ c °sI v Ljh - v £mp + <“ +,)l 
- c °s( v LJh — Vgjrip + ( w + 01 

H-sin l v Lj h - v cmp + (w + n! 


S can be found from ( 1 7) by substituting u> + f + jt/ 2 for + f where it explicitly appears in the 

above equation. (The u + f hidden in V^ mp is not to be tampered with.) Hence S will not be 
written here. 

W is more easily found than R. By ( 14) and (1 5) it is 


7rR 2 C R F«r° / a s \2 “ * ^ ~ ^ 

— ( 1 2 2 2 1j 2 d L A N01 


w = _ 

2cM 


N=0 L=0 J=0 e=0 m=0 


(2N+1) 1/2 [i-(-l) N+L+C _ 2 +L +C 

2 2 ^LJh^S^Emp^ 

<2-6ojX2-4 0 m^L+lM28 + l) n+5 0 J>< l+ W h =0 p=0 
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Wix««») ^ -,pL 

\2fi-l / \2fi+l 


N L fi-l\/N L fi-1 


0 0 0 / \0 J m-1 


/ Afi+i\ As \ / N L C+l\ / N L fi+1 

4fi-m+l )( fi -m+2) ( ) - pi ) 

\2fi+3 7 \2C+l/ l /l 

J \o 0 0 / \0 J m -1 


-sin (V Uh + vL mp _n) 

I 

+cos( v uh + v g'mp- n > 

II 

*“ (V Uh + V «»p- B> 

III 

+sin<vS uh +v Lp- Si) 

IV 


+ < V(C-m-l)(fi-m) 


N L C — 1\ / N L fi-1 


0 0 0 / \0 J m+1 

N L C+l\ /N L e+i' 


- V( £+m+l)(C+m+2) 


-sin vV^ h + V g mp + n ) 
+cos(V^j h + Vg 

mp + n ^ 
+co * v Uh + V fimp + a ) 
+sin ( v LJh + V C mp + a) 


. / b fi - 1\ / b 

+(-l) m x/(fi-m— 1 )(fi -m) [ ] -p 

\2C-1/ \2C+ly 


0 0 0 / \0 J m+1 



'N L C-l\ / N L C — 1 
^0 0 0 / \0 J -m-1 


- \/( £+ m+l)(®+m+ 2 ) (— - P 

\ 2fi +3/ \2C+1 


N L fi+l\ /N L fi+1 
0 0 0 / \0 J -m-1 
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r 


S L 

-sin <v Uh - V C mp _ n ) I 

• • “ C0S<V Uh - v «mp - n ) II 

+cos < v Uh- v Cmp - n ) III 

S L 

— sin ( V LJh — v £ mp ~ S2) IV 


+(-l) m 



£+m-l)(C+m) 



-\4 8— m+1 )(8— m+2) 



- P 


r 

S L 

~ sin (V Uh~ v 8mp + n > 
-cos(VLj h -V L Cmp + S2 ) 

+c°s ( v LJh ~ v Smp + n ) 

_sin(V LJh ' V 8mp + 

\ 


~ P 



/n L 8-l\ /n L 8-1 


\0 0 


0; 


\0 J -m+1 


>b 8 \ 


28+1 





/N L 8+1 


\0 J -m+1/ J 



I 

II 

III 

IV 


+2 




m)(8 +m) 



L 

0 


C — 1\ / N 


o /\ 0 


L 

J 



+\^8 — m+ 1 )(8 +m+ 1 ) 





L 

0 



L 

J 



+c °s( v Uh + V iLp) 

+sin ( V LJ h + v Cmp ) 
c L 

+sin(V ^ + V ) 

IJh C m P ; 


-C05( v LJh + Vp m n ) 


y 

8 mp 


I 

II 

III 

IV 




Again it is to be reminded that ( 1 7) and ( 1 8) apply only to zonal variations in albedo, i.e., K = 0 in 
these equations. 

This completes the lengthy derivation of the albedo acceleration. It remains to find the dj_ and 
bj>. 



then H has the Legendre polynomial expansion 

OO 

H = 2h l p l (cosv s ) 

L=0 

where 

H L = KIPl-1 (0)-P L+ i (0)1 

with 


P_! (0) = + 1 
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(Longman, 1962, pp. 846-847). Defining now the “daylight function” 


D= HcosJ/g= 2 cos Pl (cos 
L=0 

insures the proper behavior of the irradiance over the whole earth: ~ cos J/g in the sunlit hemi- 
sphere, ~0 in the shadowed. Using 


P L (cos\|/ s ) 


(L+l) P L+1 (cos i// s ) + LP l _j(cos^ s ) 
(2L+1) 


(Hobson, 1955, p. 32) allows writing 

OO 

D = 2 d^ Pl (cos J/«j) 

L=0 

as required, where 



and P_ 2 ( 0 ) = 0. 

The dj^ can now be found from standard tables for the Legendre polynomials (e.g., Hobson, 
1955, p. 17). A listing of the d^ for 0 < L < 6 is given in Table IV. Note that all of the d^ for 
odd L are zero, except for d j. A graph of D using the first seven terms in the summation is shown 
in Figure 4. Notice that it is an excellent approximation to the exact function. 

Finding the bg is harder work. Using the binomial expansion (e.g., Selby, 1973, p. 470) 

(1 - 2p cos t? + p)“2 = (l+q) — 2 = 1 _ 2q + 3q~ — 4q^ + . . . 

(where q = -2 p cos r\ + p^) and multiplying cos tj - p gives the following approximation B' for 

B: 
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B' = (cos 17 — p) (1 - 2q + 3q 2 - . . . ) 

«• (_p+2p 3 -3p 5 + 4p 7 -5 p 9 + 6p U -7p 13 +8p 15 -9p 17 + 10p 19 -llp 21 + 12p 23 

- 1 3p 23 + 1 4 p 27 - 1 5p 29 + 1 6 p 3 ^ ) 

+ (1 - 6p 2 + 15p 4 - 28p 6 + 45p 8 - 66p 10 + 91p 12 - 120p 14 + 153p 16 - 190p 18 + 231p 20 

- 276p 22 + 325p 24 - 378 p 26 + 435p 28 - 496p 30 ) cos t? 

+ (4p-24p 3 + 72p 5 - i60p 7 + 300p 9 - 504p n +784p 13 - 1152p 15 + 1620p 17 - 2200p 19 
+ 2904p 21 - 3744 p 23 + 4732p 25 - 588 Qp 27 + 7200p 29 ) cos 2 rj 

+ (12p 2 - 80p 4 + 280p 6 -720p 8 + 1540p 10 - 29 1 2 p 1 2 +5040p 1 4 -81 60p 16 + 12540p 8 
-18480p 20 + 26312p 22 - 36400P 24 +49140p 26 - 64960p 28 ) cos 3 r? 

+ (32 p3 _ 240p 5 + 960p 7 - 2800p 9 + 6720P 1 1 - 141 12p 13 + 26880p 15 - 47520o 17 + 7920Qp 19 

- 125840p 21 + 1 92 1 92p 23 - 283920 p 25 + 407680p 27 ) cos 4 r? 

+ (80p 4 - 672p 6 + 3024p 8 - 9856p 10 + 26208 p 12 - 60480p 14 + 125664p 16 - 240768p 18 
+432432p 20 - 736736p 22 + 1201200p 24 - 1 886976 p 26 ) cos 5 v 

+ (192 p 5 - I792p 7 + 8960p 9 - 32256p 1 1 +94C80p 13 - 236544p 15 + 532224p 17 

- 1098240p 19 + 21 141 12p 21 - 3843840p 23 + 6662656p 25 ) cos 6 r? 

4<448 P 6 _ 46 08p 8 + 25344p 10 - 99840P 12 + 316800P 14 - 861696p 16 +2086656p 18 -4612608p 20 
+ 9472320 p 22 - 1 8304000 p 24 ) cos 7 v 

+ (1024 p 7 - 1 1520p 9 + 69120P 1 1 - 295680P 13 + 1013760p 15 - 2965248p 17 + 7687680p 19 

- 18120960p 21 + 39536640p 23 ) cos 8 rj 

+ (2304 p 8 - 28 160 p 10 + 1 83040 p 2 - 844800p 14 + 31 1 1680 p 1 6 - 9737728p 18 + 26906880p 20 

- 67358720p 22 ) cos 9 rj 

+ (5 1 20 p 9 - 67584P 1 1 + 473088P 1 3 - 23429 1 2p 1 5 •: 92252 1 6p 1 7 - 30750720P 1 9 
+ 90202 1 1 2p 21 ) cos 10 rj 

+ d 1264 p 10 - 1 59744 p 12 + 1 198080p 14 - 63365 1 2p 1 6 + 26557440p 18 - 93929472p 20 )cos 1 ! ij 
+ (24576p U - 372736 p 13 + 2981888p 15 - 16773120p 17 + 74547200p 19 ) cos 12 r? 
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+ (53248p 1 2 - 860 1 60p 1 4 + 73 1 1 360p 1 6 - 4358 1 440 p 1 8 ) cos 1 3 r? 


+ (1 14688 pi 3 - 1966080 p 15 <• 17694720p 17 ) cos 14 r? 

+ (245760p * 4 - 4456448p*^) cos^ t? 

+ (524288p 15 )cos 16 t? O 9 ) 

after terminating the expansion at the term (inclusive). The approximation given above holds 
only for r? < Arc cos p ; elsewhere B' = 0. 

The bg can be found from the inner product (e.g., Hobson, 1955, p. 42): 


Hence, 


b 9 = 


bo = 


28+1 

2 

26+1 


-)/: 
■)/.’ ■■ 


B Pg(cosr?) sinr? dr? . 


Pg(cosr?) d (cost? ) 


( 20 ) 


where bg is the approximation to bg. The lower limit of this last integral reflects the fact that B' = 0 
in the region invisible to the satellite. The bg’ will be used ir. place of bg in what follows. 

It is easier to compute the bg* numerically than to give analytical expressions for them. Using 
the basic equation 


r = 


a( 1 - e“) 


1 + e cos f 

and the numerical values for R^, a, and e for Lageos result in 


( 21 ) 


p = 0.520049 ( 1 + e cos f) 


so that 

p n 2 : (0.520049) n (l + ne cos 0 

since e is small. This equation, together with (19) and (20), yield the bg listed in Table V for 6 
running from 0 to 6. The values for bg' in turn allow 
6 

B' ^ ^ bg Pg (cos r?) 

6=0 

to be computed. This function is shown in Figure 5 plotted against r? together with the exact ana- 
lytical expression for B. Both curves assume Lageos to be in a circular orbit. The figure shows 
that the approximation for B is not as good as the one for D. 
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4. Earth Albedo 

Now that the formalism for the albedo acceleration has been devcloped.it is time to consider what 
model to use for the earth’s albedo for short wavelengths. Since only the zonal terms in the albedo 
are retained (K = 0), it comes down to the matter of choosing Aqqj, ^ior ^201* etc - Stephens 
et ah (1981) summarized several years’ worth of satellite data from which a set of coefficients can 
be derived. They give the monthly zonal mean albedos for latitude bands 10 degrees wide, plus the 
yearly average for each band (see their Table 4b)- Assuming that the albedo values refer to the mid- 
point of each band (so that the value for 60°N - 70°N is >aken to be at 65°N, for example), one 
can do a least-squares fit to the data with a sum of Legendre polynomials. 

The resulting coefficients Aqqj through A401 for each month, which come from fitting the first 
five polynomials, are shown in Figure 6. This figure shows that the coefficients for even N are all 
fairly constant throughout the year. The coefficients for odd N show a distinct sinusoidal behavior 
with a period of one year. However, A 3 qj is not sinusoidal about zero like A j q j . Rather, it shows 
a constant offset which is presumably due to the presence of the large Antarctic ice cap. The coef- 
ficient Aqqj gives a rough measure of Imw good the fit is. since its time-averaged value is actually 
0.30, varying only between 0.29 and 0.32 during the year (Stephens et aL, 1981. p. 9744). The 
least-squares estimate is consistently high, but only a little. 

Clearly it will do no great violence to the data to take 

F 0 qIcos 0) 

Pj q( c °s g) 

P 2 ,o( c °s0 ) 

P 3O (cos0) 

F 4 q(cos 0 ) ( 22 ) 

as the model for the earth’s albedo. This is a hybrid expression: Aqqj is Stephens et. ah’s (1981) 
time-averaged value, while the other coefficients are based on Figure 6. Obviously a more accurate 
expression could be obtained by analyzing Stephens j£l aL’s (1981) original data, but (22) is good 
enough for the purposes considered here (see Figures 6 and 7). 

Note that the coefficients for even N are assumed to be constant. The time-varying parts of 
A j q j and A 3 qj are modeled with cos f§; this can be done due to the fortuitous *ircumstance of 
perihelion passage occurring near the first of the year. Note also that the albedo is dominated by 
the zeroth and second degree harmonics. 


A(0, X ) =0.30 

+ 10.05 cos f^l 
+ 0.10 

+ [-0.01 +0.03 cosf s ] 
+ 0.04 
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5. Lageos 

The analytical formalism and the albedo model will now be applied to the Lageos satellite. Lageos 
has the folio*. ?! parameters in addition to those shown in Table I: 

n 5: n 0 + rt (t - t Q ) 

« 0 + (t - 1 0 ) 

where 

D 0 = +28.5596 deg 
u 0 = +171.9271 deg 
t Q = 2902.5 days (Julian Date) • 

For further information about Lageos, see Smith and Dunn ( 1 980) and Rubincam (1982). 

The Lageos values will be used to compute R, S, and W via (17), (18) and the modification of 
R which yields S. Other vaiues necessary for the computation are those given in Table II, the d^ 
given in Table IV, the b^' given in Table V.andtheA^Qi 8» v en by (22 ). Also, the earth’s orbit about 
the sun will be taken to be circular, so that (a^/r^) = 1. Further, 

i^S + = + ^o + ( w + f) (t - t Q ) 

where 

(w s + f s ) 0 = 42.0087 deg. 

The subscripts will run between the following values in the summations: 0 < N ^ 4, 0 < L < 4, 
and 0 < C < 5. 

What does one get for R, S. and W using all of this data? The answer is: thousands of terms for 
each of the accelerations, for each value of N. Many of the terms are zero, since 1 — ( — I )N-+L-+C an( j 
the 3-j symbols are often zero. A computer program is necessary to make all of the calculations. 
It also proves convenient to gather together all of the terms with the same frequency with the aid 
of the computer. 

Sample calculations for R, S, and W as they vary over one revolution arc shown in Figures 8-12 
foi the major terms in the albedo N = 0 and N - 2. In all of these figures the sun is on the equator 
and Lageos passes directly over the srbsolar point. The correct qualitative behavior for R when 
N - 0 (Levin. 1962; see also Anselmo et al., 1983) is evident in Figure 8. Here the acceleration is 
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roughly zero in the shadow, rises to a maximum when Lageos is directly over the subsolar point 
(where the reflection is strongest), and then falls off symmetrically towards zero. Likewise Fig- 
ure 10 show, the correct qualitative antisymmetric behavior for S (Levin, 1962; Anselmo et a]., 
1983): nearly zero in the shadow, and then negative as the satellite approaches the subsolar point 
and the sunny area pushes Lageos against the direction of motion. S is zero because of symmetry 
when Lageos is over the subsolar point; then positive as the sunny area pushes Lageos along its way. 
Note that the maximum amplitude for S is about an order-of-magnitude less than that for R. also in 
agreement with Levin (1962). 

One way in which Figure 8 differs from Levin’s figure is that his R has a fiat top. Also, the 
humps in Levin’s S are separated by a stretch where S is zero, while the humps of Figure 10 are 
joined together, fhesc differences are due to Levin’s oversimplification of the problem. 

No figure shows W for N = 0; it is zero over the entire orbit, by symmetry. However, W for 
N = 2 is shown in Figure 1 2. 

The long-period perturbations in the Keplerian elements a, e, u, 1 and ft for N = 0 are found 
next using the Gaussian form of Lagrange’s planetary equations. Perturbations in the mean anom- 
aly M are not considered. The assumption of a circular orbit will be dropped when considering 
changes in a, e, and w but retained when considering changes in I and ft. 

The change in the semimajor axis with time is given by (c.g., Blanco ind McCuskey. 1961, p. 178) 


da 

dt 


( (er sin f)R + a( 1 - e^ )S 1 

nr^l-e 2 


( (e sin OR + ( 1 + e cos f)S 1 

nVl -e^ 


where n is the mean motion and use has been made of (2 1 ). 
Substituting R 0 ' foi R and 


S' s S Q * + (e cos f) S,' 


(23) 


in the above equation yields 


da 

dt 


2e 


l.p. n\/l -e~ 


| (sin DRo+ cos f (S^ 


+ S{ )) 


(24) 
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to first order in e for the long-period (l.p.) part of da/dt. The prime denotes the parts of R and S 
which contribute to the long-period perturbations, while R Q ' and S Q ' are the portions of R' and S' 
respectively which arc independent of the eccentricity. S' is written to first order in e; hence the 
term (c cos 0 Sj' in (23). The factor c cos f ultimately comes from the eccentric part of p. A term 

2 S 0 '/(n Vi-c^) wnich appears after making the substitution of (23) into (24) has been dropped, 
s : n<e it contributes only short-period turns. 

The long-period terms in R and S are all of those which contain ±(to+f) inside the trigonometric 
argument; all other terms are short-period. According to the computer program they are 

R q ' = 10 -12 I + 1 1.729 cos[(u>s + fij) - (to+f) ) 

- 1 1.729 cos({u >5 + f«j) + (to+f) ) 

♦ 19.829 cos((w s + f s ) - (<o+f)-fil 
+ 40.293 cos(tu) 5 + f§) + ( w + 0 -O 1 
+ 1 .735 cos[(cj 5 + fg) - (to+f)+I2) 

+ 0.854 cos( (to <j + fg) + (w+f)+fi ]} ms -2 (25) 

S Q ' = I0 ~ 12 | - ! .215 sin((u><; + fg) - (o+f) I 

- 1.215 smKcjg + fg) + (to+f) 1 

- 2.05 4 sin|(u? + .'$) - (to+f) -Hi 

H 4. 174 sin|(to g + f<j) + (w+l)-fi I 

- 0.180 sin | (to 5 + fjj) - (to+f) + Hi 

+ 0.088 sin i (to s + f s ) + (to+f) + Hl} ms ~ 2 (26) 

S j ' = I0 “* 2 | + 0.754 sin ((to $ + fg) - (to+f) 1 

+ 0.754 sin|(tj£ + fj) + (o>+0 1 

+ 1.275 sin|(to ^ + f<j) - (to+f)-fi) 

- 2.591 sinlfcj s + f s ) + (to+f) -Hi 

+ 0.1 12 sinlko s + f s ) - (to+f)^ nl 

- 0.055sin|io s + f s ) + (o.»+f)+nl} ms -2 (27) 

Substituting (25) - (27) in (24) and using (16) yields for the change in semimujor axis 

£a - Aa = - 0.40 cos(( to<j + f^) - to 1 

- 0.62cos((to<j + f<j) + to 1 

- 0.95 cosl(to<j + f<j) - to -III 

+ 3.82 cos((u>,j + f^) + to — ST I 

- 0.05 cos|(< + fg) - to + r** 

+ 0.03 cos((to<j + fg) + to + Hi J mm (28) 
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after integrating and dropping some more short-period terms. Here Aa Q is Aa at time t Q . The per- 
tMrbations in the semimajor axis are thus quite small. The largest term has an amplitude of 3.82mm 
and a period of 833 days. A graph of (28) is shown in Figure 3. 

Th«. change in eccentricity with time is given by (e.g. Blanco and McCuskey, 1961, p. 178) 


de \/l-e 2 

f (e sin O R + 

a(l-e 2 ) r 

Q 

— 


j 

dt nae 

l 

r a _ 



so that 


de 

dt 




(sin f )R 0 ' + (2 cos 0 S Q ' 


after again making use of (21 ). Substituting (25) and (26) in this equation and integrating gives 


CAe-aeJa , 0 -« 

e 


— 0.844 cosKug + fg) - cj 1 

— 1.304 cosUoj + fg) + cj 1 

— 2.000 cos((c'g + fg) — cj— n 1 

+ 8.028 cos[(wg + fg) + cj — ft 1 

0.097 cos[(w 5 + fg) — cc+fi ] 
+ 0.066 cos[(cjg + fg) + cj + £2 1 


(29) 


Thus e changes by only about one part in 10^. A graph of (29) is shown in Figure 13 also. 
Proceeding in similar fashion, 


dcj VT? 

dt - nae 


(cos OR — sin f 


[h 


r 

S + 

ersin fco+0 cot I 


*> 

all e-). 


a( 1 — e“) 


W 


(e.g. Blanco and McCuskey, 1961, p. 178), so that 


dcj 


v/ 1_e2 

[ (ccs 0R o ’ - (2 sin f) S ' J 

dt 

l.p. 

nae 

f o j 
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and the term containing W has been ignored, since it is multiplied by e. Using (25) md (26) and 
integrating yields 

Acj - Au q = — 0. 174 cos[(u>g + fg) - <jj ] 

+ 0.269 cos[(ujg I fg) + uj ] 

— 0.412 cos[(wg + fg) — cj — SI ) 

— 1.656 cos[(tjg + fg) + cj —SI] 

— 0.020 cos[(cjg + fg) — cj + fi] 

— 0.014 cos[(og + fg) + cj + fi 1 arc sec • 

See Figure 14 for? graph. 

The equations above give all of the long-period terms for the stipulated conditions on N, L, and 
C. Unfortunately, the corresponding equations for the perturbations in the inclination and longi- 
tude of node contain a great many more terms. This is became W. which controls these perturba- 
tions, contains a large number of contributing terms, unlike K .nd S. Only those terms in Al and 
An which have an amplitude greater than 10 -4 arc sec will be listed here: 

Al — Al Q 10 -4 | 1.795 sin[ 

- 0.501 sin [ 2 Zl ] 

+ 0.395 sin (2 (cj g + fg) — fi j 

+ 3.018 sin[2(cjg + f s ) - 2n 1 

+ 0.205 sin[2(ujg + fg) - 3f2 ] 

- 0.141 sin[2(cj g + fg) - 4n 1 

- 0.191 sin 14 (cj g + fg) _ 4H 1 } arc sec 

An-An o = l0 _4 { 3.956 cosi • n) 

+ 0.360 cos[ 2n 1 

- 0.257 cos[2(cjg + fg) ) 

- 0.870 cosi2(cj g + fg) - n 1 

+ 1.234 cost 2(wg + fg) - 2n 1 

- 0. 1 15 cos[2(cjg + fg) — 3n 1 J arc sec 

111 esc perturbations are small. 

These expressions are again derived from the classical equations (e.g., Blanco and McCuskey, 
1 96 1 , p. 1 78). W itself is not given. For graphs of the above perturbations, see Figures 1 3 and 1 4. 
There is also a secular contribution to the node: 
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dn A I 

= - 7.856 x 10 * arc sec yr 1 

dt sec 

which is also small. This completes the derivation of the perturbations in the Ktplerian elements 
for N = 0. 

The perturbations due to the higher harmonics (N > 0) in the albedo (22) are generally smaller 
than those for N = 0. For instance, the largest term in R Q ' for N = 2 yields an amplitude of 0.29n.m 
for the perturbation in the semimajor axis. Since the perturbations for N = 0 are already small, the 
perturbations for N > 0 will not be computed, with one exception. It turns out that the N = 2 and 

N = 4 harmonics contribute significantly to the secular motion of the node, being +8.390 x 10“^ 
and +2.3 10 x 10 - ^ arc sec yr - *, respectively. The second degree term more than cancels the zeroth 
degree term, giving a total of 

s + 2.844 x 10 - ^ arc sec yr - * 

dt | sec 

for all three terms. 

6. Discussion 

The formalism developed here results in some rather lengthy equations for the albedo accelera- 
tions. However, a virtue of these equations is that they hold whether the satellite sees the terminator 
or not; in other words, they are perfectly valid for any sun-satellite geometry. Hence there is no 
need for sepaiate equations for the terminator-visible and the terminator-invisib’e cases, as in past 
treatments. 

Application of the formalism to Lageos’ orbit shows the effect of sunlight diffusely reflected 
from a uniform earth to be quite small The periodic perturbations in the node n, for instance, are 
far below the residuals level of about 0.03 arc sec of the satellite laser ranging data (Rubincam, 
1984, Figure 2). Also, the semimajor axis a is perturbed at the 5mm level over a period of 833 
days. This is smaller than the effect of some individual solar eclipses on the orbit of Lageos. The 
eclipse of 28 'arch 1979, for example, increased the semimajor axis by about 17 6mm within the 
space of an hour (Rubincam and Weiss, 1985). The albedo perturbations in e, I, and cu are likewise 
tiny. (Interestingly, the perturbations in a, e, and are dominated by the radial acceleration R.) The 
perturbations due to the higher harmonics (N > 0) in the albedo are also small. Hence diffuse re- 
flection plays only a minor role in long-period perturbations in Lageos' orbit. 
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It is clear from the periodic nature of (28) that l^ambertian reflection cannot explain the ob- 
served secular decay of 1.1 mm day - * in the semimajor axis (e.g., Rubincam, 1982). This has been 
attributed to charged particle drag (Afonso et al , 1980; Mignard, 1981; Rubincam, 1982; A- 
fonso et al., 1985). Morgan (1984), however, suggests the decay may be caused by forward scat- 
tering due to a morning-evening asymmetry in the albedo. 

Lambertian reflection cannot explain the observed unmodeled along-track variations in accel- 
eration, either. These long-period variations have magnitude ~ 2 x 1C - *"* ms - ^ (e.g., Rubincam, 
1985; Anselmo et ah, 1983). It is obvious from { .7) that the along-track albedo acceleration is of 
magnitude e x 2.5 x 10 ms ~ = 0.01 x 10 “ ms . This is two orders of magnitude too 
small to explain the observations. 

The physical reason why the long-period variations in S depend on the eccentricity e is as fol- 
lows. Consider the effect of a single surface dement on a circular orbit. The surface element re- 
flects according to Lambert’s law. On t!..- arc where the satellite approaches the element radiation 
pressure pushes the satellite against .he direction of motion, giving a negative contribution to S. 
Correspondingly on the arc where the satellite moves away from the element the satellite is pushed 
in the direction of motion, giving a positive contribution to S. The positive and negative contribu- 
tions will exactly cancel when summed over a circular orbit because of the symmetry of the diffuse 
reflection required by Lambert’s law. The same will obviously hold true for a sum over surface ele- 
ments. Hence there can be no long-penod variations in along-track acceleration for a circular orbit; 
only short-period. An eccentric orbit must be invoked to give the asymmetry needed to produce 
long-period variations in S; thus the appearance of e. 

This argument will hold regardless of the changes in albedo across the face of the earth. There- 
fore a north-south asymmetry in the albedo cannot be invoked to explain the observed variations 
in along-track acceleration as by Anselmo et al. (1983), at least not with Lambertian reflection. 
Walch (1982) also finds Lambert's law inadequate to explain the observations. He numerically 
modeled a combination of Lambertian reflection, specular reflection, and infrared radiation using a 
zonal albedo. His resulting acceleration is about an order-of-magnitude smaller than observed. This 
is still an order-of-magnitude larger than that computed here. The difference is presumably due to 
his inclusion of specular reflection, which is not considered here. 

What then is causing the observed variations in along-track acceleration? One possibility is 
fluctuations in charged particle drag ( Afonso et al., 1985; Rubincam, 1982). Another is varia- 
tions in forward scattering (Morgan, 1984) 

How well do the assumptions made here mimic reality? There have been some obvious simpli- 
fications: higher harmonics and tesseral harmonics in the albedo, a w.. length-dependent albedo, 
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an oblate earth, a ragged terminator, etc. have all been neglected. Perhaps more serious is the as- 
sumption of a time-averaged albedo, rather than a time-varying albedo. Also, more Legendre poly- 
nomials could certainly be used to improve the approximation for B, i.e., bring the dashed line 
closer to the solid line in Figure 5. 

All of the above probably produced negligible error compared to the main assumption: that the 
earth reflects sunlight according to Lambert’s law. In fact the ways in which the earth reflects 
light are still poorly understood (Stephens et al., 1981, p. 9741). Clouds and land behave some- 
what like Lambert’s law (Jacobowitz, 1981; Stowe and Taylor, 1981; Cube, 1982, pp. 57-58), but 
the oceans display a distinct anisotropic reflectance pattern with a strong specular component 
(Stowe and Taylor, 1981). Snow and ice are also non-Lambertian reflectors (Stephens et al. 1981, 
p. 9741). Hence the main future endeavor in this field should be in pursuing laws of reflectance 
more realistic than Lambert’s law. However, the problem will become much more complicated be- 
cause of angular dependence as well as the space and time variations of reflecting surf ices. A start 
has been made in this direction with specular reflection (Walch, 1982) and forv/ard scattering 
(Morgan, 1984). 

As for the present formalism, it can be modified to examine the perturbations in Lageos’ 
orbit due to infrared radiation. These can then be compared to Sehnal’s (1981) results. This work 
is now in progress. 
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APPENDIX I 


This appendix shows that the force in Lambert’s law cannot be written as the gradient of a 
potential . 

Suppose the satellite’s position is at rj in Figure 3; the force is directed along rj and is propor- 
tional to 



where ri” = Xi^ + yi^ + Zi^. *f this force were derivable from a potential, then one could write 
•• 1 1 2 2 2 2 

r| =VV where (3V/3x) = XjZj/(xj + yj + Zj ) , etc. apart from constant multiplicative fac- 
tors. Then it should certainly be true that 


Working this out implies that 


av 

dzi 


. 3 |i- V -l 

3y, 


yi 


2\2 


= 0. 


( x j + yj + z i ) 


Since yj is not necessarily equal to zero, the above equation cannot always hold. Hence there can 
be no Lambert law potential. 

It is easy to see physically why this is so: the work done over a closed circuit is not zero. Con- 
sider taking a particle around the elementary circuit ABCD in Figure 3. The circuit is assumed for 
simplicity to lie in the y j - Zj plane. No work is done over circular arcs BC and DA since the force 
is perpendicular to the path. Over AB the work done is Arj cos apart from constant fac- 

tors where is the angle between AB and the Zj axis. Likewise the work done over CD is 
-Arj cos 02^1^’ w h ere 07 ls ^e angle between CD and the Zj axis. The net work is thus 
A r j (cos 0j - cos 0->)/r|^. This is not zero, since (Jj 4 02- Hence the force is not conservative and 
therefore not derivable from a potential. 
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APPENDIX 2 


The quantum-mechanical spherical harmonics Y L J (9,X) are related to the geophysical spherical 
harmonics Yjjj( 0, A) by the equation 


i (— 1) “ p 

Y l J (0,a) = ft 

1^(2-6 0 j)J'/2 L 


_ Lj J|1 ^ ’ x >- - 2«j|j|)i Y L)J | 2 (0.X) 


The inverse relations are 


y LJ 1^> x > = f«(2-6 oJ )j' /2 [(-l) J Y l J (0X) + Y l ~ j (0X)J , 

Iff ( 2 - 6 ,)1 1/2 

Y LJ2 (fl ’ X) = [(-1 ) J Y, J (0,X) - Y r - J (0,»| , 

i 

The geophysical spherical harmonics are defined by 
Y LJ1^’ X ) = Pl,j(cos 0 ) cos mA 
y U 2(0> x ) = ^L,J^ cos sin mA 


where 


PL,j(cos0) = (2 -5 oJ )(2L+l) — J)? ‘ /2 p 

(L+J)! J U 


(cos0) 


so that 


n TT 

Y Uj( 0 ^ ) Y L'J'j’ 


J >0 


I I Y Uj (fl ’ X ^ Y L'J'j' (0lX) sin0(10 dx = 4ff6 LL ' Sjj' 6 -'. 

J 0 *'0 JJ 

In other words, the geophysical spherical harmonics follow Kaula's (1966, p. 7)4, normalization 
ote that Y L02 («.») = 0. Likewise any coefficient C L ' of Y L Vx) is related to the coefficients 
L L|J|1- C L|J|2 of the respective Y L|J| ,(0,A), Y L|J|2 (0.A) by 
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J+IJI 

C L J = [n(2 -6 oJ )l 1/2 (-1) 2 [C lijh +(1 -26j|j 1 )iC L ,j| 2 ] 


and 



C LJ2 “ 


1 


|4„<2 -« oJ)J W 


i 


[4?r(2 - ^jj)l 


Vi 


[(-u J c L J + c L - J ] , 
[ ( -i) J c L J -c L - J ] , 


J 2*0 


J > 0 


Obviously, C|_q 2 = 0- 
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TABLE CAPTIONS 


Table I. Notation for quantities relating to Lageos. The quantities in the first part of the table re- 
fer to Lag'-os’ orbit, while those ir. the second part refer to the satellite itself. Dash (-) indicates 
various numerical values arc assumed. 

Table II. Notations for quantities relating to the sun. earth, and universe. Dash (--) indicates 
various numerical values are assumed. Symbols not shown in Tables I and II are defined in the 
text. 

Table III. Unnormalized inclinations functions Fg m p(0 for 8 = 0, 1. 

Table IV. Values for d^. 

Table V. Values for bjJ for Lageos, based on P = 0.52004 1 ) (1 * ecosf). 
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TABLE I 


Notation for quantities relating to Lageos. The quantities in the first 
part of the table refer to Lageos’ orbit, while those in the second part 
refer to the satellite itself. Dash ( - ) indicates various numerical 
values are assumed. 


Quantity 

Symbol 

Numerical Value 

Scmimajor axis 

a 

1.227 x 10 7 m 

Eccentricity 

e 

0.004 

True anomaly 

f 

- 

Inclination 

I 

1 09. 9 deg 

Mean motion 

n 

4.645 x lO^ 4 s" 1 

Earth-Lagcos distance 

r L 

- 

Argument of perigee 

CJ 

- 

Perigee rate 

U) 

0.21 12 deg day" 

Longitude of node 

n 

i 

Node rate 

n 

0.3425 deg day 

R V r L 

p 

- 

Radiation coefficient 

C R 

1.13 

Mass 

m l 

407 kg 

Radius 

r l 

0.3 m 


TABLE 11 

Notation for quantities relating to the sun. earth and universe. Dash 
(-) indicates various numerical values arc assumed. Symbols not 
shown in Tables 1 and II are defined in the text. 

Quantity Symbol Numerical Value 

Semimajor axis 

a S 

1.496 x 10 1 1 m 

True anomaly 

f S 

- 

Inclination 

*S 

23.44 deg 

Earth-sun distance 

r S 

- 

Argument of perigee 

UJg 

- 

Longitude of nude 

i “ S . 

0 

Daily motion 

‘*S ♦ fs 

0.9856 deg day I 

Solar constant 

F s° 

1376 W m" 2 

Earth albedo 

A(0.A) 

- 

Radius of albedo surface 

r a 

6381 km 

Speed of light 

c 

2 .998 x 10 B ms 1 

Gravitation constant 

G 

6.67 x 10 -11 m 3 kg 1 

Time 

t 
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TABI.F III 


Unnormali/cd inclination functions F^. (I) hr 
C=0. I. 


« 

0 

1 

1 

1 


m P PCmp ^ 0 

0 0 1 

0 0 (sinI)/2 

0 1 (-sinIl/2 

1 0 (1 + co»I)/2 

1 1 (1 -co»I)/2 


TABLE IV 
Values for d^ 


0 +1/4 

1 +1/2 

2 +5/16 

3 0 

4 -3 m 


5 0 

6 +13/256 
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TABLE V 


Values for bg' for Lageos, based on p = 0.520049 (1 + e cos f). 


t 



b e 

i 



0 

0 

+0.45795 

+ 

0.37014 e cos f 

i 

+1 .23068 

+ 

1.21735 e cos f 

2 

+1.64127 

T 

2.25323 e cos f 

3 

+1.63083 

+ 

3.20260 e cos f 

4 

+ 1 31099 

+ 

3.74560 e cos f 

5 

+0.88601 

+ 

2.91399 e cos f 

6 

+0.50452 

+ 

2.71441 e cos f 
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FIGURE CAPTIONS 


Figure l. The Lageos satellite. Lageos is a completely passive satellite whose aluminum surface is 
studded with laser retroreflectors. 

Figure 2. Positions of the sun, satellite, and surface element in the inertial coordinate system xyz. 
The coordinate system does not rotate around with the earth. The x-axis points to the vernal equi- 
nox and the z-axis lies along the rotation axis. Since the system is considered to be inertial, fic- 
titious forces produced by the earth’s motion around the sun are ignored. 

Figure 3. The local frame XjyjZj of the surface element dS. The vector r] and the circuit ABCD 
are discussed in Appendix 1 ; otherwise they are to be ignored. 

Figure 4. The function D (dashed line), for L summed between 0 and F, and the exact function 
cos \p (solid line), both plotted against i p§. D is basically the normalised solar irradiance. 

Figure 5. The approximation for the function B (dashed line), for £ summed 0 and 6 , and the 
exact function (solid line). Both are plotted against 17 . A circular orbit is assumed. 

Figure 6. The coefficients Aqqj through A^Qjfor each month. The error bars indicate one stan- 
dard deviation. Arrows indicate error bars too big to fit on the graph. The larger errors are due 
to sparser data because of shadowing at either pole. 

Figure 7. The time-averaged zonal albedo of Stephens et al. (1981) (solid line), and the time-aver- 
aged approximation used here (dashed fine). 

Figure 8 . The radial acceleration R over one revolution for N = 0. The sun is on the equator and 
Lageos passes over the subpolar point. 

Figure 9. The radial acceleration R over one revolution for N = 2. The sun is on the equator and 
I jgcos passes over the subpolar point. 

Figure 10. The along track acceleration S ov C i one revolution for N = 0. The sun is on the equator 
and Lageos passes over the subpolar point. 
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Figure 1 1. The along-track acceleration S over one revolution for N = 2. The sun is on the equator 
and Lageos passes over the subpolar point. 

Figure 12. The accele;atior. W (nonnal to the orbital plane) over one revolution for N = 2. the 
sun is on the equator and Lageos passes over the subpolar point. 

Figure 13. The top graph shows the long-period changes in Lageos semimajor axis and eccentricity. 
The changes in a and e have virtually the same shape; hence only one graph is used for both. The 
bottom graph shows changes in inclination. Both graphs are for an earth with a uniform albedo 
(N = 0)and A 00 , =0.3. 

Figure 14. The long-period changes in the longitude of the Lageos node (top) and the argument 
of perigee (bottom) with time. The secular change in f2 is omitted. Both graphs are for N = 0, 
Aqqi = 0-3. 
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z 



Figure 2. Positions of the sun, satellite, and surface element in the inertial coordinate system xyz. The coordinate 
system does not rotate around with the earth. The x-axis points to the vernal equinox and the z-axis lies 
along the rotation axis. Since the system is considered to be inertial, fictitious forces produced by the 
earth's motion around the sun are ignored. 


57 



1 



NORMALIZED SOLAR IRRADIANCE 
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D is basically the normalized solar irradiance. 



FUNCTION B FOR LAGEOS 
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plotted against tj. A circular orbit is assumed. 




JFMAMJJASOND 

Figure 6. The coefficients Aqqj through j for each month. The error bars indicate one standard deviation. 

Arrows indicate error bars too big to lit on the graph. The larger errors are due to sparser data, because 
of shadowing at eithei pole. 
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F.gu re 7. The lime-averaged ,„„al albedo of Srephena e, al. ,198 1 ) (aolid hoe,, and ,he lime-averaged approxima.ion naed here (daahed line) 
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Figure 9. The radial acceleration R over one revolution for N = 2. The sun is on the equator and Lageos passes over the subpolar point. 
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Figure 10. The along-track acceleration S over one revolution for N = 0. The sun is on the equator and Lagcos passes ovei the subpolar point. 
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Figure 1 1 The along-track acceleration S over one revolution for N = 2. The sun is on the equator and Lasers pares* over the subpolar point. 
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Figure 12. The acceleration W (normal to ihe orbital plane) over one revolution for N - 2. The sun is on the equator and Lageos passes over 
the subpolar point. 
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LAGEOS 

CHANGES IN SEMIMAJOR AXIS, 

ECCENTRICITY, AND INCLINATION 
Aooi = 0.3 

Figure 13. The top graph shows the long-period changes in the Lageos semimajor axis and eccentricity with time. 

The changes in a and e have virtually the same shape, hence only one graph is used for both. The bot- 
tom graph shows changes in inclination. Both graphs are for an earth with a uniform albedo (N = 0) 
and Aqq. = 0.3. 
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Figure 14. The long-period changes in the longitude of the Lageos node (toj. nd the argument of perigee (bottom) 
with time. The secular change in n is omitted. Both graphs are f 0, Aqq j =0.3. 
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